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Abstract 


This  is  the  Fiscal  Year  2010  annual  report  by  the  PI  to  the  US  Office  of  Naval  Research 
(ONR)  for  Grant  N00014-06- 1-0830:  “Measurements  and  analysis  of  reverberation,  target 
echo,  and  clutter.”  The  main  accomplishments  for  the  year  were:  (i)  completion  of  a  re¬ 
port  on  monostatic  reverberation  calculations  over  a  sloping  bottom;  (ii)  development  of  a 
bistatic  range  dependent  clutter  model  that  includes  towed  array  beam  patterns;  (iii)  calcu¬ 
lations  of  time  spreading  of  scattering  features  and  target  echo  and  their  effect  on  signal 
excess;  (iv)  additional  work  on  fathometer  returns  and  reverberation  at  short  times.  The 
most  significant  accomplishment  was  the  bistatic  range-dependent  clutter  model  based  on 
adiabatic  modes;  it  is  computationally  efficient  and  is  being  implemented  in  the  DRDC 
Sonar  Test  Bed  for  potential  use  by  the  Canadian  Forces. 


Resume 


Le  present  document  est  le  rapport  pour  Texercice  2010  presente  par  le  chercheur  princi¬ 
pal  a  T  Office  of  Naval  Research  (ONR)  des  Etats-Unis  pour  la  subvention  N00014-06-1- 
0830  :  ^Measurements  and  analysis  of  reverberation,  target  echo,  and  clutter  »  (Mesures 
et  analyse  de  la  reverberation,  des  echos  de  cible  et  du  fouillis  d’echo).  Les  principales 
realisations  pour  Tannee  ont  ete  les  suivantes  :  (i)  achevement  d’un  rapport  sur  les  cal- 
culs  de  reverberation  monostatique  sur  un  fond  en  pente ;  (ii)  mise  au  point  d’un  modele 
a  fouillis  d’ echos  dependant  de  la  portee  bistatique,  incluant  des  modeles  de  faisceaux 
remorques ;  (iii)  calculs  de  la  repartition  temporelle  des  caracteristiques  de  diffusion  et 
echos  de  cible  et  leurs  effets  sur  I’exces  de  signaux ;  (iv)  travaux  additionnels  sur  les  son- 
deurs  a  ultrasons  et  la  reverberation  de  courte  duree.  La  realisation  la  plus  importante  a 
ete  le  developpement  d’un  modele  de  fouillis  d’echos  base  sur  des  modeles  adiabatiques 
et  dependant  de  la  portee  bistatique.  Ce  modele  est  efficace  sur  le  plan  des  calculs  et  une 
version  initiale  pent  etre  utilisee  sur  une  plateforme  Windows  ;  une  version  Unix  est  mise 
en  cEuvre  dans  le  banc  d’essai  des  systemes  de  RDDC  et  pourrait  etre  utilisee  par  les  Forces 
canadiennes. 
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Executive  summary 


Measurements  and  analysis  of  reverberation,  target 
echo,  and  clutter:  FY10  Annual  Report  for  Office  of 
Naval  Research  Grant  N0001 4-06-1 -0830 

Dale  D.  Ellis;  DRDC  Atlantic  ECR  2010-208;  Defence  R&D  Canada  -  Atlantic; 
December  201 0. 

Background:  The  Principal  Investigator  (PI)  spends  three  months  per  year  at  the  Applied 
Research  Laboratory  of  The  Pennsylvania  State  University,  funded  by  the  US  Office  of 
Naval  Research  (ONR).  Each  year  a  short  Annual  Report  needs  to  be  submitted  to  ONR. 

Principal  results:  This  is  the  Fiscal  Year  2010  annual  report  by  the  PI  to  the  US  Office 
of  Naval  Research  (ONR)  for  Grant  NOOO 14-06- 1-0830:  “Measurements  and  analysis  of 
reverberation,  target  echo,  and  clutter.”  The  main  accomplishments  for  the  year  were:  (i) 
completion  of  a  report  on  monostatic  reverberation  calculations  over  a  sloping  bottom;  (ii) 
development  of  a  bistatic  range  dependent  clutter  model  that  includes  towed  array  beam 
patterns;  (iii)  calculations  of  time  spreading  of  scattering  features  and  target  echo  and  their 
effect  on  signal  excess;  (iv)  additional  work  on  fathometer  returns  and  reverberation  at  short 
times.  The  most  significant  accomplishment  was  the  bistatic  range-dependent  clutter  model 
based  on  adiabatic  modes.  It  is  computationally  efficient  and  an  initial  version  is  working 
on  a  Windows  platform;  a  Unix  version  is  being  implemented  in  the  DRDC  System  Test 
Bed  for  potential  use  by  the  Canadian  Forces. 

Significance  of  results:  Many  clutter  and  target  echo  measurements  have  been  made  in  the 
US/Canada/NURC  joint  experiments  “Clutter  07”  and  “Clutter  09”  on  the  Malta  Plateau. 
A  model  to  simulate  some  of  the  results  was  needed.  The  recently-developed  clutter  model 
will  enable  some  quantitative  comparisons  to  be  made  in  a  range-dependent  environment. 
The  model  is  computationally  efficient,  and  allows  direct  comparisons  with  towed  array 
beam  time  series.  It  thus  has  potential  use  for  further  development  as  a  tactical  tool  for 
operational  exercises. 

The  ONR  Reverberation  Modeling  Workshops  in  2006  and  2008  have  stimulated  work  by 
a  number  of  researchers.  In  previous  years  the  Pi’s  normal  mode  reverberation  model  was 
shown  to  be  in  good  agreement  with  predictions  from  other  models  at  intermediate  and 
long  ranges/times.  More  recently  a  ray  model  has  been  developed  to  investigate  behaviour 
at  short  times.  In  collaboration  with  other  researchers  these  short  and  long  range  effects 
are  being  incorporated  into  a  “benchmark”  against  which  other  models  can  be  evaluated. 

The  next  step,  modelling  the  target  echo,  is  critical  to  sonar  performance  prediction.  The 
Pi’s  normal  mode  model  includes  the  time  spreading,  and  initial  results  indicate  that  time 
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spreading  causes  drop  in  signal  excess  which  is  of  the  order  of  5  dB.  This  is  consistent  with 
calculations  from  other  models  presented  at  the  2010  Workshop  on  Validation  of  Sonar 
Performance  Assessment  Tools,  and  indicates  that  sonar  models  will  need  to  incorporate 
the  effects  of  time  spreading. 

Future  work:  The  “benchmark”  paper  on  one  of  the  Reverberation  Modeling  Workshop 
problems  is  almost  ready  for  submission  to  a  journal.  The  time  spreading  work  should  also 
be  incorporated  into  a  journal  article.  The  initial  version  of  the  Clutter  Model  is  ready  to 
be  implemented  into  the  DRDC  Sonar  Test  Bed;  considerable  refinement  of  the  computa¬ 
tional  engine  is  needed,  but  the  architecture  should  allow  easy  upgrades.  Model  validation 
will  focus  on  range-dependence,  target  echo,  and  benchmarks  for  sonar  performance  as¬ 
sessment. 


IV 


DRDC  Atlantic  ECR  2010-208 


Sommaire 


Measurements  and  analysis  of  reverberation,  target 
echo,  and  clutter:  FY10  Annual  Report  for  Office  of 
Naval  Research  Grant  N0001 4-06-1 -0830 

Dale  D.  Ellis ;  DRDC  Atlantic  ECR  2010-208 ;  R  &  D  pour  la  defense  Canada  - 
Atlantique;  decembre  2010. 

Introduction  :  Le  chercheur  principal  (CP)  passe  trois  mois  par  annee  a  1’ Office  of  Naval 
Research  of  The  Pennsylvania  State  University,  finance  par  le  US  Office  of  Naval  Research 
(ONR).  A  chaque  annee,  un  court  rapport  annuel  doit  etre  presente  a  I’ONR. 

Resultats  :  Le  present  document  est  le  rapport  pour  Pexercice  2010  presente  par  le  cher¬ 
cheur  principal  a  1’ Office  of  Naval  Research  (ONR)  des  Etats-Unis  pour  la  subvention 
N00014-06-1-0830  :  ^Measurements  and  analysis  of  reverberation,  target  echo,  and  clut¬ 
ter  >  (Mesures  et  analyse  de  la  reverberation,  des  echos  de  cible  et  du  fouillis  d’echo).  Les 
principales  realisations  pour  Pannee  ont  ete  les  suivantes  :  (i)  achevement  d’un  rapport  sur 
les  calculs  de  reverberation  monostatique  sur  un  fond  en  pente;  (ii)  mise  au  point  d’un 
modele  a  fouillis  d’echos  dependant  de  la  portee  bistatique,  incluant  des  modeles  de  fais- 
ceaux  remorques  ;  (iii)  calculs  de  la  repartition  temporelle  des  caracteristiques  de  diffusion 
et  echos  de  cible  et  leurs  effets  sur  I’exces  de  signaux;  (iv)  travaux  additionnels  sur  les 
sondeurs  a  ultrasons  et  la  reverberation  de  courte  duree.  La  realisation  la  plus  importante 
a  ete  le  developpement  d’un  modele  de  fouillis  d’echos  base  sur  des  modeles  adiabatiques 
et  dependant  de  la  portee  bistatique.  Ce  modele  est  efficace  sur  le  plan  des  calculs  et  une 
version  initiale  pent  etre  utilisee  sur  une  plateforme  Windows  ;  une  version  Unix  est  mise 
en  oeuvre  dans  le  banc  d’essai  des  systemes  de  RDDC  et  pourrait  etre  utilisee  par  les  Forces 
canadiennes. 

Portee  :  Beaucoup  de  mesures  de  fouillis  d’echos  et  d’echos  de  cible  ont  ete  effectuees  lors 
des  experiences  conjointes  Etats-Unis/Canada/NURC  appelees  <  Clutter  07  »  et  <  Clut¬ 
ter  09  >  sur  le  plateau  de  Malte.  Un  modele  visant  a  simuler  certains  des  resultats  etait 
requis.  Le  modele  de  fouillis  d’echos  mis  au  point  recemment  permettra  d’effectuer  cer- 
taines  comparaisons  quantitatives  pour  un  milieu  dependant  de  la  portee.  Le  modele  est  ef¬ 
ficace  sur  le  plan  des  calculs  et  permet  d’effectuer  des  comparaisons  directes  avec  les  series 
temporelles  des  faisceaux  remorques.  II  pourrait  done  etre  utilise  pour  les  developpements 
ulterieurs  a  titre  d’outil  tactique  pour  les  exercices  operationnels. 

Les  ateliers  de  I’ONR  sur  la  modelisation  de  la  reverberation  en  2006  et  en  2008  ont  stimule 
un  certain  nombre  de  chercheurs  a  realiser  des  travaux.  Au  cours  des  annees  precedentes, 
le  modele  de  reverberation  en  mode  normal  du  chercheur  principal  s’ est  avere  en  bonne 
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concordance  avec  les  previsions  d’autres  modeles  pour  les  portees/durees  de  duree  in- 
termediaire  et  de  longue  duree.  Plus  reeemment,  un  modele  de  rayon  a  ete  mis  au  point 
pour  examiner  le  eomportement  dans  les  eourtes  durees.  En  eollaboration  avee  d’autres 
ehereheurs,  ees  effets  sur  les  eourtes  et  longues  portees  sont  integres  dans  un  <  points  de 
referenee  »  avee  lequel  d’autres  modeles  peuvent  etre  e values. 

La  proehaine  etape,  qui  eonsiste  a  modeliser  I’eeho  de  la  eible,  est  eritique  pour  la  prevision 
du  rendement  d’un  sonar.  Le  modele  normal  du  ehereheur  prineipal  eomprend  la  repartition 
temporelle,  et  les  resultats  initiaux  montrent  que  la  repartition  temporelle  entrarne  une 
baisse  des  exees  de  signaux,  de  I’ordre  de  ~5  dB.  Cela  est  eonforme  aux  ealeuls  effeetues 
a  I’aide  d’autres  modeles  qui  ont  ete  presentes  lors  de  1’ atelier  de  2010  sur  la  Validation  des 
outils  d’ evaluation  du  rendement  des  sonars,  et  nous  permet  de  eonelure  que  les  modeles 
de  sonars  devront  integrer  les  effets  de  la  repartition  temporelle. 

Recherches  futures  :  Le  doeument  sur  les  <  repere  »  portant  sur  I’un  des  problemes  sou- 
leves  lors  de  1’ atelier  sur  la  modelisation  de  la  reverberation  est  presque  pret  a  etre  presente 
a  un  journal.  Les  travaux  sur  la  repartition  temporelle  devraient  egalement  etre  integres  a 
un  artiele  de  journal.  La  version  initiale  du  modele  de  fouillis  d’eehos  est  prete  a  etre  mise 
en  oeuvre  dans  le  bane  d’essai  des  sonars  de  RDDC ;  un  peaufinement  eonsiderable  du  mo- 
teur  de  ealeul  est  requis,  mais  I’arehiteeture  devrait  permettre  des  mises  a  niveau  faeiles. 
La  validation  du  modele  portera  essentiellement  sur  la  dependanee  a  la  distanee,  sur  I’eeho 
de  la  eible  et  sur  les  reperes  destines  a  revaluation  du  rendement  des  sonars. 
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1  Long-Term  Goals 


The  long-term  goal  of  this  work  is  to  better  understand  and  model  reverberation  and  elut- 
ter  in  shallow  water  environments,  and  to  develop  teehniques  for  Rapid  Environmental 
Assessment  (REA)  and  environmentally  adaptive  sonar. 

2  Objectives 


The  eurrent  projeet  is  a  joint  eollaboration  between  Defenee  Researeh  &  Development 
Canada  -  Atlantie  (DRDC  Atlantie)  and  the  Applied  Researeh  Eaboratory  of  The  Penn¬ 
sylvania  State  University  (ARE/PSU)  to  analyze  and  model  reverberation,  target  eeho, 
and  elutter  data  in  shallow  water.  It  allows  the  Prineipal  Investigator  (PI)  to  spend  ap¬ 
proximately  three  months  eaeh  year  at  ARE/PSU.  The  eollaboration  leverages  programs 
in  Canada,  US,  and  joint  researeh  projeets  with  the  NATO  Undersea  Researeh  Centre 
(NURC).  The  primary  effort  is  analysis  and  interpretation  of  data,  together  with  devel¬ 
opment  and  validation  of  improved  modeling  algorithms. 

3  Approach 


The  PI  spends  three  months  per  year  at  ARE/PSU,  eondueting  joint  researeh  primarily  with 
Drs.  John  Preston  and  Charles  Holland.  Additional  eollaboration  takes  plaee  throughout 
the  year.  The  main  objeetive  of  this  eollaboration  is  to  analyze,  model,  and  interpret  data 
reeeived  on  towed  arrays  during  reverberation  and  elutter  sea  trials.  The  primary  outputs 
of  the  eollaboration  are  manuseripts  for  joint  publieations  in  eonferenee  proeeedings  and 
refereed  journals.  Seeondary  outputs  are  improved  models  and  algorithms. 

Eoei  of  this  eollaboration  are  Joint  Researeh  Projeets  (JRPs)  between  NURC,  Canada,  and 
several  US  researeh  laboratories  (ARE  in  partieular).  The  present  JRP  “Charaeterizing 
and  Redueing  Clutter  for  Broadband  Aetive  Sonar”  is  essentially  eomplete.  A  new  JRP 
“Modeling  and  Stimulation  for  ASW  Aetive  Sonar  Trainers”  for  years  2011-2013  is  just 
starting  up,  with  the  PI  as  the  overall  Point  of  Contaet. 

Over  the  past  few  years,  the  ONR  Reverberation  Modeling  Workshops  [1-3]  have  been  a 
focus  for  collaboration.  The  PI  extended  and  exercised  two  of  his  models  on  a  number 
of  problems  [4],  and  collaborated  with  Preston  in  developing  a  Matlab-based  model  [5]. 
The  model  is  also  being  validated  against  more  computationally-intensive  physics-based 
models  developed  by  other  researchers.  Other  collaborations  with  Holland  and  Ainslie  led 
to  presentations  at  a  Special  Session  of  the  Acoustical  Society  of  America  in  November 
2007  [6, 7].  A  paper  with  Holland  has  been  published  [8],  and  the  paper  with  Ainslie  is  in 
progress.  Sub-bottom  scattering  from  the  collaboration  with  Holland  [8]  is  now  an  option 


DRDC  Atlantic  ECR  2010-208 


1 


in  the  NOGRP  model.  The  PI  was  a  member  of  the  problem  definition  eommittee  for  the 
seeond  Reverberation  Modeling  Workshop  held  in  May  2008. 

Reeent  work  has  foeused  on  range-dependent  reverberation  modeing  and  target  eeho  eal- 
eulations,  both  of  whieh  are  important  for  developing  a  “Clutter  Model”. 

4  Work  Completed 


The  ONR  Reverberation  Modeling  Workshops  have  stimulated  further  work  along  the  same 
lines.  One  sueh  aetivity  was  the  2010  Symposium  on  Validation  of  Sonar  Performanee 
Assessment  Tools,  sponsored  by  the  UK  Institute  of  Aeousties  in  memory  of  David  Weston. 
One  of  the  “Weston  Symposium”  problems  extended  the  ONR  problems  to  the  full  sonar 
seenario,  ineluding  matehed  filter  proeessing,  baekground  noise,  and  signal-to-noise  ratio. 
The  main  organizer  was  M.  Ainslie  of  TNO  in  the  Netherlands,  and  the  PI  provided  adviee 
on  the  problem  definitions,  whieh  are  now  being  published  formally  [9]. 

During  the  past  year  the  PI  worked  on  the  target  eeho  problem  from  the  2008  ONR  Re¬ 
verberation  Modeling  Workshop  and  the  2010  Weston  Symposium.  While  multipath  time 
spreading  is  not  very  important  for  diffuse  reverberation,  it  ean  have  have  a  signifieant 
effeet  for  feature  seattering  and  target  eeho  [10, 11].  A  presentation  was  made  at  the  Oe- 
tober  2009  meeting  of  the  Aeoustieal  Soeiety  of  Ameriea  [12],  and  at  the  Weston  Sympo¬ 
sium  [13].  An  example  is  shown  in  the  “Results”  seetion  below. 

The  ONR  Workshops  uneovered  differenees  between  model  predietions  at  short  times;  e.g., 
ranges  less  than  about  10  water  depths,  where  the  fathometer  returns  and  steep-angle  rever¬ 
beration  are  important.  A  problem  that  generated  a  great  deal  of  diseussion  was  Workshop 
Problem  VI,  for  the  rough  bottom  with  the  summer  sound  speed  profile.  The  Pi’s  normal¬ 
mode  results  were  extended  to  short  times  and  submitted  to  the  organizers  for  eomparison 
with  other  mode  models,  and  ray-model  predietions.  They  were  ineluded  in  a  presentation 
at  the  Oetober  2009  meeting  of  the  Aeoustieal  Soeiety  of  Ameriea  [14]. 

It  is  diffieult  to  properly  inelude  the  short  range  effeets  in  a  normal  mode  model,  so, 
to  investigate  the  short-time  reverberation,  a  refinement  of  a  1986  straight-line  ray  traee 
model  [15]  was  developed  and  implemented.  This  also  required  a  eareful  look  at  the  pulse 
definition  for  the  Reverberation  Modeling  Workshop.  The  results  were  submitted  to  the 
ONR  Workshop  organizers,  and  presented  at  the  Weston  Symposium  [16]. 

As  part  of  the  model  development  for  model-data  eomparisons,  the  fast  shallow-water  re¬ 
verberation  model  [4, 10]  based  on  normal  modes  has  been  extended  to  a  bistatie,  range- 
dependent  elutter  model  that  ineludes  target  eeho  and  feature  seattering  [17-19].  The  for¬ 
mulation  based  on  adiabatie  modes  had  been  developed  earlier  [17]  and  ealeulations  were 
done  using  Matlab  for  a  flat  bottom  environment.  In  2010  the  model  was  eonverted  to 
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Fortran  95,  and  extended  to  handle  range-dependent  bathymetry  and  scattering.  Like  the 
reverberation  model,  it  is  computationally  efficient  and  includes  the  3-D  effects  of  towed 
array  beam  patterns  in  order  to  facilitate  comparison  with  experimental  measurements.  A 
short  paper  will  appear  in  Canadian  Acoustics  [19].  Some  results  are  illustrated  in  the 
“Results”  section  below. 

A  very  useful  feature  of  the  clutter  model  is  that  it  is  a  subroutine,  rather  than  a  standalone 
program.  This  facilitates  its  inclusion  in  other  applications.  Work  is  progressing  to  include 
it  in  the  DRDC  System  Test  Bed  (a  version  of  which  goes  to  sea  with  the  Canadian  Forces); 
a  standalone  version  with  Java  GUI  and  public  domain  databases  is  already  running  on  a 
Windows  system  [20].  Further  model  improvements  will  occur  in  the  next  year,  apply¬ 
ing  it  to  other  problems  from  the  ONR  Workshops  and  Weston  Symposiums,  as  well  as 
comparisons  with  towed  array  data,  particularly  from  the  Clutter  ’07  and  Clutter  ’09  trials. 


5  Results 


This  section  illustrates  a  few  examples  from  activities  during  the  past  year. 

Time  spreading  of  scattering  features  and  target  echo 

Figure  1  illustrates  the  reverberation  from  a  flat  bottom  with  Lambert  scattering.  The  envi¬ 
ronment  is  as  described  in  Problem  IX  of  the  Reverberation  Workshop.  The  problem  was 
modified  by  increasing  the  scattering  strength  by  10  dB  from  -27  dB  to  -17  dB,  for  a  100  m 
annulus  between  20.0  to  20.1  km.  If  the  travel  time  for  all  the  paths  is  assumed  to  be  the 
same,  as  is  typically  done  in  many  formulations,  the  reverberation  would  jump  by  10  dB 
for  0.133  s.  The  blue  curve  shows  the  computations  done  using  a  normal  mode  reverber¬ 
ation  model  at  250  Hz.  If,  however,  the  different  multipaths  take  different  travel  times, 
then  the  scattering  from  the  feature  is  spread  in  time;  the  green  curve  shows  the  impulse 
response  using  the  modal  group  velocities,  and  the  red  curve  shows  the  impulse  convolved 
with  a  uniform  pulse  of  duration  0.08  s.  (Note  that  the  red  curve  overlies  much  of  the 
green  curve.)  Outside  the  scattering  region,  the  curves  with  and  without  the  group  velocity 
correction  differ  by  about  0.5  dB.  This  is  because  the  constant  group  velocity  was  chosen 
to  be  1500  m/s,  but  the  actual  group  velocities  are  all  less  than  this;  thus,  the  scattering  at  a 
given  instant  of  time  is  coming  from  shorter  ranges  where  the  propagation  loss  is  less,  and 
so  gives  slightly  higher  reverberation.  The  small  difference  is  generally  not  important,  so 
the  extra  computation  to  handle  the  different  modal  travel  times  is  not  usually  needed  for 
reverberation  calculations. 

However,  the  situation  is  much  different  for  the  target  echo,  as  illustrated  in  Fig.  2.  If  all 
the  multipaths  arrived  at  the  same  time,  then  the  echo  would  have  the  same  shape  as  the 
pulse,  and  the  peak  would  lie  on  the  dark  green  curve.  The  echo-to-reverberation  ratio  for 
the  example  shown  would  be  about  -3  dB.  However,  with  different  multipath  travel  times 


DRDC  Atlantic  ECR  2010-208 


3 


10  dB  enhanced  scattering  from  20.0  to  20.1  km  (0.133s) 


Figure  1:  Reverberation  structure  at  250  Hz  from  a  scattering  annulus  enhanced  by  10  dB 
for  100  m  extent  in  range. 


Prob  A2I,  250  Hz,  Range  10  km,  tau0=0. 0941 38s 


Figure  2:  Echo  time  structure  for  a  point  target  at  range  1 0  km. 
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Echo-to- Reverberation;  src-rec.  separation  20  km 


X  (km) 


Figure  3:  Echo-to-reverberation  ratio  in  100  by  100  km  area  using  ideal  beam  patterns. 
The  source  is  indicated  by  a  *,  and  the  receiver  by  a  o  with  heading  225°  indicated  by  the 
short  line;  the  time  ellipses  indicate  20,  40  and  60  s. 

the  echo  is  stretched  out  in  time,  and  the  amplitude  is  reduced.  The  light  green  dots  in  the 
figure  correspond  to  the  different  arrival  times  of  the  mode  pairs,  and  the  red  curve  is  the 
Hann-shaped  pulse  convolved  with  the  impulse  responses.  In  this  case  the  peak  response  is 
reduced  by  about  6  dB,  so  the  peak  echo  to  reverberation  ratio  now  about  -9  dB.  This  can 
be  quite  important  for  detection  predictions  based  on  signal  excess. 

These  results  shown  are  for  normal  modes,  but  similar  results  can  be  obtained  with  other 
models  that  include  time  spreading  (such  as  Holland’s  energy  flux  approach). 

Range  dependent  clutter  model 

The  formulation  for  a  range-dependent  clutter  model  based  on  adiabatic  modes  had  been 
developed  earlier  [17]  and  calculations  done  for  a  flat  bottom.  In  2010  the  initial  version  of 
the  bistatic  range-dependent  version  became  operational.  Two  calculations  are  illustrated 
in  Figs.  3  and  4;  these  calculations  will  appear  in  a  short  paper  in  Canadian  Acoustics  [19]. 

Figure  3  shows  the  signal  excess  (target  echo-to-reverberation  level  in  dB)  on  a  grid  gen¬ 
eralized  from  Fig.  2  in  [17].  The  area  is  100  km  square  with  coordinates  between  {x,y)  = 
(-50,0)  km  and  {x,y)  =  (50,100)  km.  A  50x50  grid  is  used  for  the  calculations,  so  the  cen¬ 
tre  of  the  grid  points  range  from  approximately  (-49,1)  km  to  (49,99)  km  with  increments 
of  2  km  in  each  direction.  The  water  depth  is  100  m,  except  for  two  ridges  in  the  y-direction 
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Time  (s) 

Figure  4:  Towed  array  reverberation  beam  time  series  corresponding  to  Fig.  3. 

rising  to  60  m,  and  another  ridge  rising  to  70  m  in  the  x-direetion;  both  have  gaps  near  the 
middle.  There  is  a  single  seamount  of  height  50  m  near  {x^y)  =  (-40,85)  km.  The  bottom 
has  Lambert  seattering  with  a  strength  of  -27  dB,  exeept  for  a  -I-IO  dB  enhancement  along 
the  line  (2,2)  km  to  (50,50)  km.  Similarly  the  target  (at  depth  10  m)  has  echo  strength  of 
8  dB,  except  for  a  7  dB  enhancement  along  the  line  (-48,52)  km  to  (2,2)  km.  The  source  is 
at  (-10,45)  km  at  depth  30  m,  and  the  receiver  at  (10,45)  km  at  depth  50  m. 

The  basic  environment  is  similar  to  the  ONR  3D  problems.  It  has  isospeed  water  of 
1500  m/s  over  a  sand  bottom  half  space  of  density  2.0  relative  to  seawater,  sound  speed 
1750  m/s  and  attenuation  of  0.5  dB/wavelength;  the  volume  absorption  in  the  water  is  a 
version  of  Thorp’s  formula.  The  source  is  omni-directional  with  unit  energy  (10  dB  source 
level  for  a  duration  of  0.1  s),  operating  at  a  frequency  of  250  Hz.  The  towed  array  was 
chosen  to  give  a  horizontal  beam  width  of  3.6°;  39  omnidirectional  elements  at  spacing  of 
2.5  m  with  Hann  weights  were  used  for  the  beam  time  series.  The  CPU  time  on  a  2  GHz 
computer  was  only  1.5  s,  with  modes  (usually  16)  being  calculated  at  each  of  the  2500  grid 
points. 

For  the  grid  calculation,  the  beam  pattern  corresponds  to  an  ideal  one  sided  broadside 
beam  steered  at  all  azimuths;  this  gives  a  nice  picture  of  the  physical  scattering.  In  Fig.  3 
the  receiver  sees  high  reverberation  (low  echo-to-reverberation  ratio)  in  the  direction  of 
the  source  (west,  W);  similarly  along  the  high  scattering  line  to  the  SE.  The  higher  echo 
strength  along  the  line  to  the  SW  shows  up  clearly  too.  Along  the  three  ridges,  the  echo 
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to  reverberation  ratio  first  drops  (due  to  the  higher  reverberation  on  the  up  slope),  then  in¬ 
creases  on  the  down  slope.  Beyond  the  ridges  one  would  expect  some  shading  due  to  mode 
cutoff,  but  none  is  observed;  perhaps  the  reverberation  and  echo  are  affected  similarly.  The 
single  point  to  the  NW  affects  the  echo-to-reverberation  ratio  in  the  adjacent  cells.  The  el¬ 
lipses  correspond  to  20  s,  40  s  and  60  s  after  the  short  pulse,  and  are  useful  for  identifying 
features  in  the  beam  time  series  shown  in  Fig.  4. 

For  the  beam  time  series,  the  actual  conical  beam  pattern  with  left-right  ambiguity  is  used, 
so  it  can  be  directly  compared  with  measured  towed  array  data.  However,  interpretation 
is  not  as  straightforward  as  for  Fig.  3.  Figure  4  shows  time  series  corresponding  to  Fig.  3, 
and  a  selection  of  beam  steering  angles  relative  to  the  towed  array  heading  of  225°.  The 
predictions  generally  seem  to  make  sense.  The  reverberation  on  the  omni  receiver  is  10- 
15  dB  above  the  beam  predictions;  at  short  times  the  reverberation  on  the  45°  and  60° 
beams  is  higher  since  they  look  in  the  direction  of  the  source  (note  the  beams  have  left-right 
symmetry);  the  45°  beam  also  seems  to  be  picking  up  backscatter  from  the  ridges  to  the 
W  (at  about  25  s),  the  ridge  to  the  S  (at  about  40  s),  and  from  the  strong  scattering  line  (to 
the  south  at  53  s);  the  0°  and  359°  beams  are  essentially  identical;  the  90°  and  270°  beams 
should  be  identical,  and  have  the  lowest  reverberation,  except  in  the  region  of  features;  at 
long  ranges  the  endfire  beams  (0°  and  180°)  should  have  the  highest  reverberation,  and  in 
a  uniform  environment  approach  each  other. 

The  model  continues  to  evolve.  Presently  the  same  sound  speed  profile  and  bottom  loss  are 
used  at  all  locations.  The  next  step  will  be  to  generalize  it.  For  production  calculations,  it 
makes  sense  to  pre-calculate  the  modes  on  the  grid,  perhaps  interpolating  them  to  a  finer 
scale.  Then,  for  other  source-receiver  geometries  and  multistatic  scenarios,  they  can  be  re¬ 
used  for  rapid  reverberation  and  target  echo  calculations.  The  model  is  being  incorporated 
in  the  DRDC  System  Test  Bed;  a  stand  alone  Windows  version  with  a  Java  interface  is 
already  working  [20]. 

6  Impact/Applications 


From  an  operational  perspective,  clutter  is  viewed  as  one  of  the  most  important  problems 
facing  active  sonar  in  shallow  water.  The  long-term  objective  of  this  work  is  to  better 
understand  and  model  reverberation  and  clutter  in  shallow  water  environments,  and  to  de¬ 
velop  techniques  for  Rapid  Environmental  Assessment  (REA)  and  environmentally  adap¬ 
tive  sonar.  Parts  of  the  earlier  research  have  been  spun  off  into  a  DRDC  TIAPS  (Towed 
Integrated  Active-Passive  Sonar)  Technology  Demonstrator  which  has  been  evaluated  in 
ASW  exercises  against  submarine  targets.  The  work  on  clutter  is  related  to  the  DRDC 
effort  in  auralization  and  co-operative  work  with  TTCP  and  other  ONR  efforts. 

If  the  target  echo  model  can  be  validated,  this  could  be  a  useful  method  for  estimating 
the  target  strength  of  clutter  features — and  even  submarines — in  multipath  shallow  water 
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environments. 


One  goal  is  to  be  able  to  use  the  model  with  real  elutter  data  from  a  towed  array.  One  eould 
subtraet  out  the  baekground  reverberation,  ineluding  range-dependent  effeets  and  known 
seattering  features,  leaving  behind  the  unidentified  elutter  on  a  display.  These  unidentified 
features  would  then  be  investigated  by  other  teehniques  to  try  to  determine  their  nature. 

7  Transitions 


The  range-dependent  reverberation  and  target  eeho  model  is  being  implemented  as  the 
Clutter  Model  on  the  DRDC  System  Test  Bed  (a  prototype  sonar  proeessing  system),  a 
version  of  whieh  is  the  Pleiades  System  used  in  some  Canadian  Forees  exereises.  Researeh 
eontraets  to  ineorporate  the  reverberation  and  target  eeho  were  let  to  General  Dynamies 
of  Canada  Limited  in  2009,  and  a  follow-on  eontraet  to  Brooke  Numerieal  Serviees  in 
2010.  A  standalone  version  with  publie  domain  databases  and  a  Java  GUI  is  presently 
operational  [20]  under  a  Windows  operating  system,  and  integration  with  the  System  Test 
Bed  is  antieipated  early  in  201 1. 

The  David  Weston  Sonar  Performanee  Assessment  Symposium  was  held  in  Cambridge, 
UK,  7-9  April  2010.  It  had  number  of  seenarios  based  on  the  ONR  Reverberation  Mod¬ 
eling  Workshop  problems,  extended  to  the  eomplete  sonar  problem.  The  driving  foree 
behind  the  sonar  modeling  is  the  Low  Frequeney  Aetive  Sonar  program  of  TNG  and  the 
Royal  Netherlands  Navy. 

8  Related  Projects 


This  projeet  has  eontributed  to  the  US/Canada/NURC  Joint  Researeh  Projeet  “Charaeteriz- 
ing  and  Redueing  Clutter  in  Broadband  Aetive  Sonar”  whieh  reeeived  substantial  funding 
from  ONR.  A  new  proposal  “Modelling  and  Stimulation  for  ASW  Aetive  Sonar  Trainers” 
has  been  approved  for  the  201 1-2013  Seientifie  Program  of  Work  at  NURC. 

This  ONR  projeet  also  eontributes  to  the  DRDC  Atlantie  researeh  program: 
http://www.atlantio.drdo-rddo.go.oa/researohteoh/soienoe_seotions_eng.html, 
in  partioular.  Underwater  Sensing, 

http://www.atlantio.drdo-rddo.go.oa/researohteoh/us_eng.html. 

As  well,  the  personal  interaotion  on  this  projeet  faoilitates  additional  oollaborations  be¬ 
tween  soientists  in  the  various  researeh  laboratories. 
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